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We study the work fluctuations of two types of finite quantum spin chains under the application 
of a time-dependent magnetic field in the context of the fluctuation relation and Jarzynski equal- 
ity. The two types of quantum chains correspond to the integrable fsing quantum chain and the 
non-integrable XX quantum chain in a longitudinal magnetic field. For several magnetic field pro- 
tocols, the quantum Crooks and Jarzynski relations are numerically tested and fulfilled. As a more 
interesting situation, we consider the forcing regime where a periodic magnetic field is applied. In 
the Ising case we give an exact solution in terms of double confluent Heun functions. We show that 
the fluctuations of the work performed by the external periodic drift are maximum at a frequency 
proportional to the amplitude of the field. In the non-integrable case, we show that depending on 
the field frequency a sharp transition is observed between a Poisson limit work distribution at high 
frequencies toward a normal work distribution at low frequencies. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

The study of fluctuations in non-equilibrium small sys- 
tems has become an active field of research during the last 
years. The reasons are twofold: on one hand, nanoscaled 
systems are nowadays quite easily manufactured, opening 
the door to the emergence of various nano-technologies. 
On the other hand in the field of non-equilibrium statis- 
tical mechanics, where exact results are very few, the dis- 
covery of fluctuation symmetries [l], expressed by the 
Gallavotti-Cohen fluctuation theorernjj, 13 [H, [|| 0, H| and 
Jarzynski equality [9( have opened new theoretical per- 
spectives. The fluctuation theorem is a statement on the 
time-reversal symmetry of the fluctuations of the entropy 
production along a non-equilibrium path Q|. Whereas 
this theorem is an asymptotic statement, Crooks derived 
an interesting identity reading [H[ 
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AS is the entropy production for a system driven during 
a time r within a forward protocol (t), Pp is the dis- 
tribution of entropy production in the forward process 
and P B is the distribution associated to the backward 
process, that is when the system is driven in a time- 
reversed manner. The celebrated Jarzynski equality @ 
is easily derived from the Crooks relation. Indeed, uti- 
lizing for a thermalised system AS" = f}(W — AF), with 
AF = [F(X f (t)) - F(A(0))] and integrating over W we 
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which is Jarzynski relation. These relations, initially de- 
rived for classical systems, have been extended to the 
quantal world as well [H, [ij, [3 [H, QJ, [l7| within var- 
ious setup and prescriptions for the actual measurement 
of work performed on a quantum system [bSj . 

As it is clear from the extensivity of the entropy pro- 
duction, the probability to observe a decrease of entropy 
in a given non-equilibrium trajectory is exponentially 
small for a macroscopic system. However, for micro- 
scopic or mesoscopic systems these untypical trajectories 
arise with a significant probability and consequently can 
be observed and measured in actual experimental setup. 
One may mention in particular experimental tests on the 
stretching of RNA molecules [HI, [2(| , experiments on tor- 
sional pendulum [2l| , on colloidal particl es J22I [23| , on 
photochromic defect center in diamonds [241 ]. See Ref. 
[2 a ] for an excellent review of the experimental related 
investigations. It is then interesting to calculate or pre- 
dict the shape of these fluctuations in some explicit mod- 
els for small system size [26|, [2?J . In this contribution 
we focus our attention to the work fluctuations of two 
time-dependent quantum spin chains where the time de- 
pendence is due to an externally applied magnetic field. 
The two models considered are the integrable quantum 
Ising chain in a time-dependent transverse field and the 
XX-quantum chain with a longitudinal magnetic field 
that breaks its integrability. We consider these two cases 
as archetypical of the integrable and non-integrable sit- 
uations and we expect that the work distributions will 
reflect somehow these differences. We check Jarzynski 
and Crooks relations and compute explicitly the non- 
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equilibrium work distributions associated to several mag- 
netic field protocols. In particular, we study the limiting 
steady distributions in the important case of periodically 
driven system. 

The paper is organized as follows: in the next section 
we define the models and protocols for studying work 
fluctuations. We insist in particular on the definition of 
work we are using in our study, pointing explicitly as 
it was done in [l|| HI] the discrepancies between this 
definition and the use of a work operator. We present 
then the numerical results corroborating Jarzynski and 
Crooks relations before turning to the main part which 
is the study of work fluctuations within a periodic drift. 
An exact solution, in terms of double confluent Heun 
functions is given in the Ising case. The study of the 
XX-cham is only numerical, but nevertheless evidences 
are drown out to show clearly the appearance of a sharp 
transition between an exponential work distribution at 
high field frequencies toward a gaussian work distribu- 
tion at lower frequencies. We summarize and discuss our 
results in the last section. 



II. QUANTUM SPIN CHAINS IN A TIME 
DEPENDENT FIELD 

A. Definition of the model and protocol 

We study quantum chains with a time-dependent field, 
in the context of fluctuations of the work performed on 
them by the time-dependent force, defined through the 
following Hamiltonian: 

N-l .N-l N 

^ = - 2 E «+i - 2 E °i °W - 2 k{t) E < ( 3 ) 

i— 1 i— 1 i— 1 

in particular in the two cases A = and A = 1. The 
<t's are Pauli's matrices and h(t) is a time-dependend 
magnetic field applied in the z-direction and leading to 
the Zeeman term —h(t)M z where M z is the total z- 
component of the magnetization. 

At A = 0, it corresponds to the Ising quantum chain. 
In a static transverse field the Ising chain can be mapped 
after a Jordan- Wigner transformation onto a fermionic 
problem which can be diagonalised after a suitable canon- 
ical (Bogoliubov) transformation. In the thermodynamic 
limit N — » oo, this model presents a quantum phase 
transition at h = h c = 1 which is in the 2d classi- 
cal Ising model universality class. The integrability and 
non-triviality (in the sense of physical properties) of the 
model is at the origin of its wide use in many fundamen- 
tal studies concerning in particular the testing grounds 
of non-equilibrium quantum statistical mechanics. One 
may mention for example relaxation properties from an 
inhomogeneous initial state [2!| G^l l3ll. I&j l33l . l34j , en- 
tropy production after a local quench [35l |36I |. entropy 
production in the unique NESS (non-equilibrium steady 
state) generated from a two temperature initial state 



[37| , or even relaxation properties after a quench through 
the critical point (basically h is varied from hi > h c to 

h } < h c ) da GUI. 

The case A = 1 corresponds to the isotropic XX- 
quantum chain in a time-dependent magnetic field, which 
breaks the integrability of the model. With a periodic de- 
pendence of the field, this model has to be compared to 
the quantum kicked Heisenberg chain: 

H ^ = -\ E + «+i + s ® Vt %<ri) ( 4 ) 

where the time-dependent perturbation is periodically 
switched on. This model was studied in the context of the 
transition from integrability to ergodicity in the thermo- 
dynamic limit [4(| IU ■ Here instead of the two-body 
interaction (in the fermionic picture) of of +1 we apply a 
magnetic field in the z-direction which leads to a many- 
body fermionic term. 

In the following, we consider the distribution of the 
work performed by the time-dependent field h(t) from 
the initial time ti to the final time tf. The varying field 
will lead to transitions between the initial state of the 
system to a new state at the final time. Here we use as a 
definition for the work the difference of energies in the fi- 
nal state and initial state, Wif = AEif = Ef —Ei, which 
assumes that we have measurements performed at U and 
tf. The work distribution P(W) is a weighted sum over 
all initial and final states of the quantum transition prob- 
ability \{4>f, Utfii ) | 2 , where U is the unitary time-evolution 
operator associated to the time-dependent Hamiltonian 
H(t), and is given by 

P(W) =J2^E if - W)\(cf>f, Uc/>i)\ 2 (5) 

where <fn and </>/ are respectively initial and final eigen- 
states of the initial and final Hamiltonian and uj^ >i — 
((pi, po4>i) where po ex e~^ H ^ is the initial (which is sup- 
posed to be canonical here) density matrix. It is clear 
that this definition differs from the introduction of a work 
operator defined such that its expectation in the state p 
would give the actual performed work: 

W = / dsTr{p(s)^—} (6) 

which can be shown to be given by 

W = Tr{p{ti)8H} = 
Tr{p(U)[U+(t f ,U)H(t f )U(t f ,U) - H(U)}} (7) 

after a short algebra. One is then tempted to define 
as a work operator 8H = H(tf) — H(ti), with A the 
Heisenberg picture of A. However, it is clear that in gen- 
eral the moments ((SH) n ) would differ from the moments 
(W n ) = JdW W n P(W) with P(W) defined in © and 
consequently the expectation of the operator e _/3<5ff 

(e-" w ) = Tr{p(t r )e-f )SH } (8) 
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will not necessarily give the Jarzynski result . Nev- 

ertheless, one can show that the moments so defined will 
differ only for n > 2, so that for the average and average 
square work one can use the work operator SH since then 

(SH) = Tr{p(U)SH} = J dW WP{W) 

{(SH) 2 ) = Tr{p(U)(SH) 2 } = J dW W 2 P(W) (9) 

with P(W) = J2ijS(AE if - W)\((f>f ,U<f>i)\ 2 u+*. The 
equality of the first moments is trivial. To prove the 
equality for the second moments, one may notices that 
for [p(ti), H(ti)] — using the cyclically property of the 
trace one has Tr{p(ti)[H(tf), H (ti)]} — 0, which then 
leads trivially to the equality of the second moments: 
JdWW 2 P(W) = Y.ij{Ef - E % ) 2 U h U&)\ 2 u** = 
Tr{p(U)(H(t f ) - H(U)) 2 } - Tr{p(U)[H(t f ),H(U)}} = 
Tr{p(ti)(H(tf) — H(ti)) 2 }. In short, the second moments 
agree whatever the commutator [H (t f) , H (ti)] is. At the 
level of the third moments, n = 3, one may show that a 
non-vanishing term Tr{p(ti)[H(tf), H(ti)]H(tf)} comes 
into play leading to a difference between the moments 
associated to the two preceding definitions of the work. 

In order to compute the work distribution one needs 
to know the unitary time-evolution operator U for the 
given field protocol, which is a difficult task in general. 
The numerical study assumes that the true dynamics is 
well approximated by the step like unitary evolution 

M-l 

U(t,t ) = 11 U(t n+1 ,t n ) (10) 

71=0 

with 

t/(t„,*„-i) = e ^(*»> At (ii) 

where At is an elementary time increment, meaning that 
the real function h(t) is approximated by a step func- 
tion. The true dynamics is obtained by taking the limit 
At — > 0, and M — » oo. In the present study we use 
typically values of At of order 10~ 2 (time is measured in 
units where h = 1 and the spin-spin coupling J x in the re- 
direction of the spin chains is set to one) . The transition 
probabilities = \(<pi, U(t, to), 4>f)\ 2 i where (j>i (0/) 

are eigenstates of the initial(final) Hamiltonian, are com- 
puted numerically using (fT0|) and (jTTJ) on small chains 
of typically less than 10 spins and with open boundary 
conditions. 



B. Jarzynski and Crooks relations 

In order to test our numerical procedure we compare 
first the free energy differences obtained from the Jarzyn- 
ski equality in a linear varying field protocol with 
the equilibrium free energy differences calculated directly 



from the partition function 

F(h(t)) = -~ \nTr{e- 0H{ ^} . (12) 




h 



FIG. 1: Free energy difference obtained from the Jarzynski 
equality compared to the exact value for the Ising chain with 
N — 7 as a function of the terminal field value at f3 = 1 and 
h(t) = t. 

In figure 1 we have plotted the free energy differences 
obtained from the Jarzynski equality and from a direct 
equilibrium calculation as a function of the final field 
value h at inverse temperature (3—1. The agreement 
is excellent with a relative deviation which is less than 
10~ 6 . 

We have also considered the ratio Pf(W)/ Pb(-W) = 
P(W) / P(—W) for the symmetric protocol h(t) — 
h a sin(7rt/r) for t € [0,r]. In figure 2 we see explicitly 
the exponential dependence e@ w for h a = 1/2 and 2 for 
a chain of size N — 8. Again, the agreement of the 
numerics with the analytical Crooks and Jarzynski re- 
lations is fulfilled with a relative deviation which is less 
than 10~ 6 . We have tested also the Crooks relation in 
the non-symmetrical protocol whose results are presented 
in figure 3. Again, the exponential dependence e^ w is 
clearly seen on the graph. 

C. Stationary distributions in the driven regime 

1. Ising quantum chain 

a. Numerical study We consider now the case where 
the system is forced with a periodic external field with 
period r, maximum value h and minimum value 0. We 
are interested in the limiting stationary work distribu- 
tion (if any) obtained after many periods of the external 
driving field. 

Before going on the driven situation, we concentrate 
on the fluctuations of the work after one period. In fig- 
ure 2] we show the distributions obtained on the Ising 



4 




W W 

FIG. 2: Crooks fluctuation ratio in the symmetrical case for FIG. 4: Work probability distribution for the Ising model at 
the XX chain at A = 1 with N = 8 at = 1 and r = 1. with N = 7 spins at /3 = for h a = 1/2. 
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FIG. 3: Crooks fluctuation ratio in the asymmetrical case for 
the XX chain at A = 1 with N = 8 at (3 = 1. 



chain with N = 7 spins for h — 1/2 and for differ- 
ent periods r of the forcing field for an initial infinite 
temperature state p = (equidistribution of the initial 
eigenstates). In that case, obviously, the average work is 
vanishing, which is reflected in the symmetry of the work 
distribution P(W) = P(—W). On figure[4]we see clearly 
that the evolution of the distribution is not monotonic 
with the forcing period r. The width of the distribution 
first increases with r, reaching a maximum and then de- 
creases with t. This behavior is clearly seen on the sec- 
ond moment (W 2 ) in figure [5] where we show the variance 
(W 2 )/3=o as a function of r for different maximum field 
values h for a N = 7 spin chain. We see that the width of 
the distribution increases with the amplitude of the per- 
turbation, since then the time-dependent perturbation is 
more effectively coupled with the spin chain. Moreover, 
the period T max associated to the maximum width de- 
creases as the amplitude is increased with, as seen from 
the numerics, a law T rnax ~ l/h for large fields enough 
(h > 1). Moreover we have seen numerically that T max 
is almost independent of the system size. As the field 



is increased, the maximum value of the variance reaches 
an asymptotic finite value which is linearly depending 
on the size N of the chain. One has from the numerics 
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FIG. 5: (W 2 ) for the Ising model with N = 7 spins at f3 = 
as a function of the period r. 

We have also computed the first and second moments 
of the distribution as a function of r for different inverse 
temperatures p. The initial state temperature acts as a 
scale factor for the moments: 

(W)/3 = tanh(/3/2)(cW0 (13) 

where (5W) , defined by this equation, is the temperature 
independent part of the average work. As seen on figure 
[5] one has for the variance 

(W 2 ) - (W) 2 ~ (W\ (14) 

meaning that the width of the distribution has a very 
small dependence on the initial state temperature, at 
least in the considered range of p. 

We turn now to the driven situation with many peri- 
ods. In figure[7]we present the time evolution of the aver- 
age work for different periods r at a field value h Q = 1/2 
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FIG. 6: Variance of the work distribution for the Ising model 
with N — 7 spins at h = 1 as a function of the period r for 
P = 0, P = 1/2 and @ = 1. 



and inverse temperature (3 = 1. As expected, for large 
periods the system follows almost adiabatically the vari- 
ation of the associated equilibrium free energy (the re- 
versible work). Moreover, it is seen as expected from the 
second principle that the actual average work is always 
greater than the free energy difference (see the inset). 
As the period of the external transverse field is lowered, 
we see that the deviation from the free energy difference 
becomes larger and one is no more able to recognize the 
underlying oscillatory external field. The largest vari- 
ations and most "chaotic" behavior is obtained at the 
value T max previously identified for one period. At higher 
frequencies the fluctuations of the average work decrease 
again. 




FIG. 7: Time evolution of the average work for the Ising 
chain with a triangular periodic transverse field at /3 — 1 and 
maximal amplitude h — 1/2 with periods t = 1 (circles), 
t — 3 (squares) and r = 5 (triangles). The inset shows the 
free energy difference (lower curve) and the average work (top 
curve) for a large period r = 50 of the field. In this case the 
transformation is nearly adiabatic. 



b. Exact solution for the driven Ising model With a 
sinusoidal forcing field h(t) oc sm(oji) it is in fact possible 
to solve exactly the dynamics of the Ising quantum chain. 
Indeed, following the lines of Ref . [43| . we can put the 
periodic boundary conditions Hamiltonian ([3]), after a 
Jordan Wigner transformation followed by a Bogoliubov 
one, into a sum of commuting operators: 



with 



N/2 

P =i 



H p = (cos0 p - h(t))[a,pa p + al p a_ p ] 
.1 . 



-t-sm 



t> p [cLp al p + cipd-p] + hit) 



(15) 

(16) 
(17) 



where [H p , H q ] = and where 



2-wp/N. The a" 1 



are Fermi operators in momentum space. Using the fol- 
lowing basis (|0), a+atp 1 0), a+ 1 0), a± p | 0}), where |0) is 
the vacuum of the a fermions, we have the 4x4 rep- 
resentation 



h(t) —ism<j)p 

ism<j)p 2cos0 p — h(t) 

coscj) n 




•p 

cos 



(18) 



In the Heisenberg picture, the time evolution matrix 
U p (t) in the subset p is governed by the equation (H = 1) 



i±Up(t) 



Up{t)H p (t) 



(19) 



with the boundary condition Z7 P (0) = /. Using (fT8| and 
(|19[) . it is easy to obtain for the non vanishing elements 
of the 4x4 matrix U p (t) the differential equation 



U?{t) U™(t) 



and (7 p 33 (i) = U^{t) 



dt v K\t) uf(t) 

h(t) 

i sinri 



-% sm< 



2 cos < 



p 

h(t) 



(20) 



From the two-coupled 



first order differential equations, it is easy to obtain the 
decoupled second order ones. For example, one has for 
the 11-component 

iu" -2 cos <pu' ~(ti -ism 2 cf) + ih[2 cos <f)-h])u = (21) 

with u(0) = 1, ii'(0) = — ih(0) and similar equations 
for the other components. Taking for the field the form 
h{t) = h\ + h>x cos(cjt) and putting it in the differential 
equation (|2"Tj) one gets solutions in terms of double con- 
fluent Heun function Hd, namely for Up X (t) and Up 2 one 
has 

1 



U l 1{ ti = [(a»^-Me M *H D (a,/3,7,*,-is) 

-h 1 )e- iA m I ,{-a,(3, 1 ,5,iz)\ ,(22) 



2h 2 

- (cos (j> p — 2hi 



r; 



-e- lAt UD{~a,P,- 1: 5,iz)] , (23) 

with z = tan(wi/2), A t = — sm(ujt) — tcos(f>p, a = 
7 /2 = 4/i 2 /", = 4, [/1 2 - |/i 2 - ftj + e^| 2 ] and 5 = 
— -% [/i 2 , — |/i2 + fci + e'^l 2 ] , and similar expressions for 
the other components. 

The equilibrium initial state is factorized as p(0) cc 
rip^i e~P^p(o) and its time-evolution is given by the ten- 
sor product 

p(t) = pi(t)p 2 (t)...p N/2 (t) (24) 

with 

p p (t) = U p (t)p p (p)U+(t). (25) 

The average work performed on the system during the 
time t with the periodic forcing is given by the integral 
of the magnetization: 

(W)(t) = - [ ds h'(s)Tr{p(s)M z } (26) 
Jo 

where M z = J^pli^p^ with M z = a+a p + at p a- p — 1, is 
the total magnetization operator in the z-direction. Us- 
ing the factorised form of p(t) and the additive structure 
of M z , it is easy to compute (M z ){t) as a sum of N/2 
independent modes. One obtains 

N/2 taml ( PK \ 

{M z )(t) = £ l^l[(/ l( o)_cos^)(2|C/ p 11 | 2 -l) 

p—1 y 

- 2sin0 p /m{t/ p 12 t/ p 11 *}] (27) 

with A p = \J sin 2 <f> p + (h(0) — cos <f>p) 2 - Together with 
the solutions of the U p components in terms of Heun 
functions, this formally solves the problem for the average 
work. One may notice that the factor tanh( ^4p^ ) in this 
expression leads to the previously observed tanh(/3/2) in 
(fT3| since in that case h(0) = and then A p (0) = 1 
and one has for the average magnetization 

(M z )(t) = - tanh (t\[l + 2j2 cos 4> P | Uf | 2 

\ / p 

+ sin0 p /m{y p 12 [/ p 11 *}] (28) 

and consequently the average work can be written as 
in (HU): {W) f3 = tanh(/3/2)(<5W^). The independence 
of the A^/2 modes also explains the linear behavior 
of the variance observed in the previous section since 
(W*) - (W)} = E" /2 (W*)p - (W q )} ~ N, where W q 
is the work associated to the gth mode. Finally, taking 
hi = —hi = h /2, such that h(t) — h sm 2 (^), we have 
an exact solution for a situation which closely resembles 
to the triangular drift treated above numerically. We 
show the behavior of the average work on figure [5] for 
different periods r = 2ix /lo. As previously observed, the 
average work fluctuations are maximum around r = 3. 
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FIG. 8: Time evolution of the average work for the Ising chain 
with a sine square periodic transverse field at j3 = 1 and 
maximal amplitude h = 1/2 with periods r = 1 (circles), 
t = 3 (squares) and r — 4 (triangles). The inset shows the 
free energy difference (lower curve) and the average work (top 
curve) for a large period r = 50 of the field. In this case the 
transformation is nearly adiabatic. 



2. XX chain in a longitudinal magnetic field 

Finally we present the results obtained in the forced 
regime for the XX chain. In figure [5] we show the sta- 
tionary work probability distributions for short and long 
period r on a chain of N — 8 spins obtained from an 
initial infinite temperature state after applying the field 
over typically several tens of periods (after about 30 pe- 
riods the work distributions are collapsing nicely on the 
same curve). We see clearly the appearance of two differ- 
ent regimes depending on the time scale r. For very short 
periods t«1, the shape of the distribution is very well 
approximated by an exponential law e~^ w ^ a /(2a) while 
for long enough time scale r its shape is close to a normal 
law. One may notice on figure [HI the very good collapse 
of the distributions for different values of the period r 
in both short and long periods regimes. We have also 
seen that this transition survives at finite temperature. 
The main difference between the infinite and the finite 
temperature cases being that in the last case the work 
distribution is no more symmetric but squeezed toward 
the positive work values as seen in figure 1101 

To give evidences of this transition, we have plotted 
on figure [TT1 the variance (W 2 ) — (W) 2 (the average work 
is vanishing for (3 = 0) of the distribution P(W) as a 
function of the time scale r. For a system of size N = 8, 
we observe a sharp increase of the variance from values 
around 0.25 at time scales smaller than for t c ~ 1 to 
values around 2.5 at larger periods r where the shape is 
gaussian-like. Note that the threshold value t c is depend- 
ing on the system size, decreasing as the size is increased. 
On figure [T2l we show the field amplitude dependence of 
the variance (W 2 )q of the work distribution for short and 
long periods of the external drift. One may notice that 



7 




-20 -10 



W 

FIG. 9: Work probability distribution for the XX quantum 
chain with N — 8 spins at j3 — for h = 1/2. 
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FIG. 10: Work probability distribution for the XX quantum 
chain with N — 8 spins at j3 = 1/2 for h = 1/2. 

at long periods, the increase of the variance with the 
field amplitude is almost linear, while it is more curved 
at short periods. On figure [T3l we see that as the size is 
increased the variance at large time scales strongly grows 
with N, while there is almost no dependence on the size 
at high frequencies. It means that when the period of 
the external field is to short, the system has no time 
to follow the perturbation and only transitions between 
nearby levels are significantly induced. For the density 
of work w = W/N this leads to a S(w) distribution in the 
limit h — > 0, since as seen on figure [12] a = yj (W 2 }/2 is 
a vanishing function of the amplitude of the field h, pre- 
sumably linear. On the contrary, and as the system size 
is increased, at long periods enough, meaning that the 
external field drives efficiently the system from its initial 
level to a new state, transitions between far apart levels 
are becoming significant. Taking that the variance seems 
to be const. hN 2 from figures [T2l and [T3l the distribution 
tt(w) in that case behaves as e ~ const - w l h . 

Finally, we have analyzed the temperature dependence 
of the work distribution P(W) in the large r regime. As 




FIG. 11: Variance of the work distribution for the XX quan- 
tum chain with N = 8 spins at f3 = and = 1 for h = 1/2 
as a function of the field period r. 
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FIG. 12: Second moment of the work distribution for the XX 
quantum chain with N — 8 spins at /3 = as a function of 
the field maximal amplitude h . 

seen on figure 1141 the average work after a first increase 
with (3 seems to finally saturate at a value close to 0.44. 
We have checked that this shape is compatible with a 
tanh(/i/3) behavior. Contrary to the almost temperature 
independent variance of the Ising case distributions, we 
see here that the variance of the distribution is decreas- 
ing(increasing) with /^(temperature), which is what one 
would normally expect. 



III. SUMMARY AND DISCUSSION 

In this study we have presented the numerical and 
analytical results we have obtained for the work dis- 
tribution of small quantum systems driven by an ex- 
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N 

FIG. 13: Second moment of the work distribution for the XX 
quantum chain with h — 1/2, at (3 = as a function of the 
system size N for r = 0.1 (squares) and r = 10 (circles). 




P 

FIG. 14: First (circles), second (squares) and variance (trian- 
gles) of the work distribution for the XX quantum chain with 
N = 8 spins at r = 10 at h = 1/2 as a function of the inverse 
temperature f3. 

ternal field. We have considered two different models, 
namely the intcgrablc Ising quantum chain in a time- 
dependent transverse field and the XX quantum chain 
in a time-dependent longitudinal magnetic field. In this 
last model, the presence of the longitudinal field breaks 
the free fermionic structure of the chain, while the trans- 
verse field in the Ising case preserves the free particles 
structure which obviously leads to integrability. In a 
first stage, we have checked the validity of both quan- 
tum Jarzynski and quantum Crooks relations that were 
largely discussed in section 2. After this initial check, 
validating our numerical approach, we focused our at- 
tention to the periodic driven situation. In the Ising 
chain, exploiting its free fermionic structure, we solved 



exactly the unitary dynamics of the system in terms of 
somehow complicated double confluent Hcun functions. 
The main features observed are, as expected, that for a 
slow enough process, that is large period of the oscilla- 
tions of the field compared to the coupling constant, we 
recover the adiabatic situation where the work slightly 
differs, from above, from the equilibrium free energy dif- 
ference. This is a statement of the second principle. As 
the process is fastened, the average work starts to deviate 
significantly from the free energy difference with fluctu- 
ations in time that are growing as the frequency of the 
field is increased. Nevertheless, at very high frequencies, 
the work fluctuations decrease again to zero. In this case, 
the variation of the field is so fast that the system is not 
able to follow it anymore. The maximum amplitude of 
the work fluctuations is obtained at an intermediate pe- 
riod which is of order of the inverse field amplitude, at 
least for large field enough. This threshold is understood 
as a dynamical resonance. In the Ising chain, one is not 
able to observe a steady work distribution. The behav- 
ior of the XX chain in a longitudinal periodic field is 
quite different. We observe numerically, that the long- 
time work distribution reaches a steady shape which is 
strongly dependent on the frequency of the field. In- 
deed, at high frequencies the distribution is well fitted 
by an exponential curve, e~\ w \l a /(2a), where a <~ h is 
almost size- independent. Lowering the frequency, there is 
a sharp transition toward a Gaussian like behavior for the 
work distribution with a variance which seems from the 
numerics to be proportional to hN 2 , . In order to under- 
stand the transition between these two limit distribution, 
one has to realize that the energy band of the chain has 
a width Se of order N (since the typical coupling J is set 
to one). The periodic perturbation introduces a typical 
energy scale lo = 2it/t. For lo > Se, there is no resonant 
coupling of the system with the periodic forcing. But as 
soon as lo < Se resonant coupling leads to the appearance 
of resonant peaks, that is sharp increases of the transition 
probability in the regions W — ±w and integer multiples 
of lo. The typical width of those peaks is of the order of 
the amplitude h of the perturbation. Consequently, the 
deviation to the initial exponential distribution starts at 
periods r > 0(2ir/N). As the period is increased further, 
the first peaks moves toward the center of the distribu- 
tion and new resonant peaks enter into play from the 
boundaries. Finally, the superposition of these resonant 
peaks for r > 0(1/ h) leads to a new limit distribution of 
the work. So the transition region between the two limit 
distributions is 0(2tt/N) < t < 0(1/ h). Finally, one 
may remark that the change of shape from exponential- 
like to Gaussian-like could be linked to the integrability 
of the model. Indeed, at very short periods, the system 
is not able to follow the external perturbation and its be- 
havior is governed by the initial integrablc Hamiltonian 
(XX-chain without field) while at larger periods the sys- 
tem feels effectively its non-integrability, possibly leading 
to a Gaussian distribution of level spacing (which is our 
work). To confirm eventually this scenario one needs to 
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push further this investigation. 
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